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section turns out to be very sensitive to the value of the saturation scale. 
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1 Introduction 

In recent years, a lot of work has been devoted to the understanding of the 
parton distributions inside highly energetic hadrons or nuclei (see for a 
pedagogical introduction). A central question in this field is the issue of parton 
saturation at small values of the momentum fraction x ||(||^]. Indeed, the 
BFKL equation (|J^,^1, when applied to extremely small values of x, leads to 
cross-sections that are too large to be compatible with unitarity constrains. 
It has been argued that this description becomes inadequate when non-linear 
effects in QCD become important, which happens when the phase-space density 
reaches the order of l/a s ^jii^. 

A theoretical description of saturation is provided by the model introduced 
by McLerran and Venugopalan [ ^2p^ , p^ | , in which small- a; gluons inside a fast 
moving nucleus are described by a classical color field, which is motivated by 
the large occupation number of the corresponding modes. The color field is 



determined by the classical Yang-Mills equation with a current induced by the 
partons at large values of x. The distribution of hard color sources inside the 
nucleus is determined by a functional density which is assumed to be Gaussian in 
the original form of the model. This physical picture has been named "Colored 
Glass Condensate". When quantum corrections are taken into account, the 
density distribution evolves according to a functional evolution equation p5| , ^6| , 
which reduces in special cases to the BFKL equation. 
It is interesting to note that this evolution equation was recently reobtained in 
the dipole picture Q , which is the natural description in the frame where the 
nucleus is at rest. 

The classical version of the color glass condensate model contains only one 
parameter, called "saturation scale" Q s , and defined as the transverse momen- 
tum scale below which saturation effects start being important. This scale 
increases with energy and with the size of the nucleus |Q0Jl3]. Therefore, at en- 
ergies high enough so that Q s 3> A QCD , the coupling constant at the saturation 
scale is small, which makes perturbative approaches feasible. 

Note that another model inspired by saturation has also been developed by 
Golec-Biernat and Wusthoff in order to fit the HERA data [^|^7j . There, sat- 
uration is included via a phenomcnological formula for the dipole cross-section 
cr(rj_) = c (l — exp(— r\/ R(x) 2 )) that saturates for large dipole sizes. The 
radius R(x) is the dipole size above which the saturation occurs, and is taken 
to be of the form R(x) — (x/xo) A ^ 2 GeV _1 . A fit to inclusive data at small 
x gave do = 23mb, xq = 3.10~ 4 and A = 0.29. This phenomenological model 
has also been successfully applied to the description of diffractivc HERA data 
[ p8[ . The main differences between this model and the colored glass condensate 
model are two-fold: in the model by Golec-Biernat and Wusthoff, saturation is 
introduced by hand, and the x evolution of the saturation scale is also ad hoc. 
In the colored glass condensate, saturation emerges explicitly as a consequence 
of solving the fully non- linear classical Yang-Mills equations, and the quantum 
evolution is described by an evolution equation that may lead in principle to 
an x dependence much more complicated than a power law. One may also note 
the attempts 29,3(1 to fill the gap between the phenomenological approach of 
Golec-Biernat and Wusthoff and perturbative QCD, by using perturbative QCD 
as a starting point and adding shadowing corrections a la Glauber-Mueller. 

Saturation effects, and in particular the color glass condensate model can be 
tested in deep inelastic scattering experiments by measuring structure functions 
like F2 f3~i| . Interestingly, it also leads to potentially observable effects in the 
field of heavy ion collisions. For central collisions, the colored glass condensate 
model has been used to compute the distribution of the gluons produced initially, 
both analytically ]3^ , |33] , p4[ and numerically on a 2-dimensional lattice |35| , |3"6| | . 
This initial gluon multiplicity is expected to control directly the observed final 
multiplicity. Recently, we have shown |37| that the parton saturation leads to 
interesting effects also in peripheral collisions, where it can be probed via the 
photoproduction of quark- antiquark pairs. In particular, the saturation scale 
Q s controls directly the position of the maximum in the transverse momentum 
spectrum of the pairs. 
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The present paper is an extension of this previous work, where we study 
the diffractive photoproduction of quark- antiquark pairs in peripheral heavy 
ion collisions. By diffractive production, we mean events in which the nuclei 
remain intact, and with a rapidity gap between the nuclei and the produced 
objects characterizing the topology of the final state. This implies that the 
exchange between the color field and the quark pair is color singlet, and carries 
a vanishing total transverse momentum. The simplest diagram for this process 
is represented in figure |l|. We consider large ultra-relativistic nuclei, and we take 




Figure 1: Prototype of the diagrams contributing to the diffractive photopro- 
duction of a qq pair in A A collisions. 

into account the electromagnetic interaction to lowest order in the coupling (or, 
equivalently, in the leading logarithmic approximation), and the interactions 
with the strong color background field to all orders. We find that the cross- 
section of this process is very sensitive to the value of the saturation scale in 
the colored glass condensate. Note that we use a Gaussian distribution in order 
to describe the hard color sources that are generating the classical color fields. 
This means that no quantum evolution has been included in this calculation. It 
is known that in the colored glass condensate model, the quantum corrections 
lead to a modification of the distribution function of hard sources when the 
energy increases p^| , p^| , p T 7|JT^ , p^ , pO| , pT|j2^ , p^ , p^ | (a similar evolution equation has 
been derived by Kovchegov and Levin in the dipole language pq ], and has 
been recently shown to be equivalent to the evolution equation in the colored 
glass condensate by Mueller ^5|). In general, this distribution will not keep a 
Gaussian form due to these corrections. However, in order to assess very simply 
saturation effects in diffractive ultra-peripheral collisions, we do not include this 
refinement here. Note that in principle, extending our calculation to include 
quantum evolution is just a matter of changing the corr elator (U(0)U f (x ± ))_ 
that appears in our results (see Eq. ([33]) for instance) and contains all the 
information about the properties of the colored glass condensate by one that 
has been calculated with the evolved distribution of "hard sources" . 

The structure of this paper is as follows. In section |2|, we summarize results 
relevant for the calculation of particle production in a classical background 
field. We also recall how to calculate diffractive quantities in the McLerran- 
Venugopalan model. In section ^, we compute the pair production amplitude 
to lowest order in the electromagnetic coupling constant. In section ^, an ex- 
pression for the transverse momentum spectrum of the components of the pair 
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is derived, and in section g results for the integrated diffractive cross-section are 
presented. Finally, sections is devoted to concluding remarks. 



2 Reduction formulas for pair production 

In this section, we remind the reader of reduction formulas one can use to com- 
pute various physical quantities in a classical background field. The formalism 
was derived in |59|, and here we are only going to recall the results, using simi- 
lar notations. Different from the situation studied in |p9fl , the relevant classical 
background field is now a superposition of an electromagnetic field and a color 
field. The color field is associated to the covariant gauge distribution of hard 
sources p (see section 2 of |$7) ) . In order to calculate physical quantities, one has 
to average over an ensemble distribution of these sources. Quantities obtained 
before averaging over the sources will be denoted by an index p. 



2.1 n versus Pi 



Certain quantities of interest for the problem of quark-antiquark production in 
a background field can be related to quark propagators with different boundary 
conditions in this background field |59) . The main results are the following. 

The probability to produce exactly one qq pair in a background field can be 
related to the Feynman (or time-ordered) propagator of a quark via 



Pi[p\ = \(0 out \O in ) m \ 



d 3 q d 3 p 
(27r) 3 2w q (27r) 3 2w p 



u(q)T®(q, -p)v(p) 



(1) 



In this equation, T F is the interacting part of the Feynman propagator of a 
quark in the background field generated by the source p. In terms of the exact 
Feynman propagator G^) and the free Feynman propagator G F , it is defined as 



Up — iap 



u p 1 p Up 



(2) 



Note in Eq. (|l]) the presence of the square of the vacuum-to- vacuum amplitude 
(0 O ut|0in)[p], which is essential for the unitarity of the resultjj]. 

Conversely, the average number of qq pairs produced in a background field, 
n — ^2nP n , where P n are the production probabilities, can be expressed in 
terms of the interacting part of the retarded propagator of a quark in this 
background: 



d 3 t 



d 3 p 



u{q)T [ ^{q, -p)v(p) 



(3) 



(2tt) 3 2w 9 (2tt) 3 2u; p 

Another difference to Eq. (0) is the absence of the vacuum-to-vacuum factor. 

'In a field theory with a time-dependent background field, the vacuum-to- vacuum ampli- 
tude is not a pure phase, but has a modulus different from one. This is a mere consequence 
of the fact that such a background field can produce or destroy particles. 
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For the background field of two nuclei in peripheral collisions, an interesting 
aspect is that in the ultra-relativistic limit the retarded propagator has a much 
simpler structure than the Feynman propagator |39| . This enables one to write 
the retarded propagator of a quark explicitly, while this is not possible for the 
time-ordered one. 



2.2 Diffractive and inclusive quantities 

So far, we have given the expression for quantities that would be observed for the 
background field configuration associated to a specific source p. However, this 
is not what is measured in experiments where quantities are usually averaged 
over many collisions. In the colored glass condensate model, this translates into 
an average over the sources, which are assumed to have a Gaussian distribution 
in the classical version of the model. 

The prescription to average observables, like for instance P\ or n, over p is 

(0)~ = J [dp\w[p\0[p\ , (4) 

with the distribution 

r~| _ / I A -J2 Pa{x~,X ± )p a (x-,X ± )\ 

Mp]=exp|-y dx dx ± 1 . (5) 

This averaging procedure yields inclusive physical quantities, "inclusive" mean- 
ing here that the final state of the nucleus which interacts by its color field is not 
observed. For the present case of qq photoproduction, it means that the process 
under consideration is: ZA —> ZqqX (the nucleus interacting electromagneti- 
cally, denoted by Z, remains intact, while the other nucleus may fragment). 

It was also noted that one can obtain a diffractive quantity (i.e. for events 
where the nucleus A does not break up) by averaging the amplitude over the 
sources p before squaring it in order to get a probability. This procedure has 
been outlined in [pZi|]4l|l and justified in Jl(|0. The difference between the two 
ways of performing the average is illustrated in figure ||, where the shaded blob 
indicates an average over the color sources. This average has the property that 
the gluons attached to one such blob carry a zero total transverse momentum 
and no net color. Given this, it is intuitive that the nucleus remains intact in the 
final state (i.e. on the cut represented by a dashed line) for the right diagram. 
In addition, this condition leads to a rapidity gap between the nucleus and the 
qq pair. More formal considerations on this issue can be found in p^] . 

Applied to the production of one qq pair, the inclusive probability is given 

by 

^-/ wHWH i 111 -"^ \-uMT<%,- pHp) \ 2 ) f . (6) 

2 For another discussion of the averaging procedure we refer the reader to [Q. 
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Figure 2: Inclusive (left) vs. diffractive (right) photoproduction of qq pairs. The 
shaded blob represents the average over a hard color source. For the diffractive 
cross-section, one has to average the amplitude before squaring it. 



Note that the vacuum-to- vacuum factor must be included inside the average, be- 
cause it depends on the source p. On the other hand, the diffractive probability 
is given by 



pdiff = 



d 3 p 



(2n) 3 2u q (2tt) 3 2cj p 



(Oo 



P m ) m u{q)T^{ q ,-p)v{p))„ 



(7) 



Regarding, on the other hand, the average number of pairs in diffractive events, 
it is important to realize that this quantity cannot be obtained from the ex- 
pression (||), which is already a combination of probabilities (i.e. squared am- 
plitudes). From Eq. (^|), therefore, only the average number of pairs produced 
in a collision, regardless of what happens to the nucleus, can be calculated J37J. 
In the rest of this paper, we focus on the calculation of P^ lS . 



3 Expression of P\ to order (Za) 2 

The results of the previous section are valid to all orders in the coupling con- 
stants occurring in the problem, in particular the electromagnetic coupling Zae q 
of a quark with charge e q (in units of the positron charge) to the nucleus 
with atomic number Z. However, the Feynman propagator and the vacuum- 
to-vacuum amplitude are not known to all orders in Za, even in the ultra- 
relativistic limit (see J39|). In this section, we show that results can be obtained 
in closed form if one works at the lowest non-trivial order in Za, but to all orders 
in the interactions with the strong color field of the other nucleus. Although 
the value of Za might seem too large (Za « 0.6 for a gold nucleus) as a useful 
expansion parameter, it turns out that at leading order the production proba- 
bilities contain large logarithms whose argument is the center of mass energy of 
the collision. It can be shown that higher order terms in Za are suppressed by 
at least one power of this large logarithm, which justifies the approximation for 
large enough collision energies. 



3.1 Pair production amplitude 

Many simplifications occur when one considers only the leading order in Za. 
Indeed, since 



u(q)n Pi (q,-p)v(p) = 0{{Za) 
|(O out |O in )~| 2 = l + 0((Za) 2 ), 



(8) 



one can approximate the vacuum-to- vacuum factor by 1 in the expression for P\ 
at leading order^ in Za. Another notable simplification occurs in the Feynman 
propagator of the quark in the background field. From an infinite series, it 
simplifies to the sum of only three terms that are represented in figure |3| The 
interaction with the color field of the second nucleus is accounted for to all 
orders by eikonal factors depicted by the black dot. Thus, one can explicitly 





> ? > 



f I 



Figure 3: The three diagrams contributing to Tf at lowest order in the elec- 
tromagnetic coupling constant. The black dot denotes the time-ordered eikonal 
matrix Tp that describes the interaction of a quark or antiquark with the colored 
glass condensate. 



write down the time-ordered scattering matrix T F in terms of the individual 
scattering matrices off the two nuclei. The sum of the three diagrams of figure 
H is formally 



T, 



QCO _ n QED 
p pi C 



QED QCD 
p Upl 17 



T 



t p — ± F '~rp ± F T -L p F F ' ± F 

where G F is the free time-ordered quark propagator 

]ji + m 



'-'F 1 F 



GIT, 



F ± F 



G%(P) 



p2 _ m 2 



(9) 



(10) 



n QCD 



is the time-ordered scattering matrix of a quark in the color field of the 



second nucleus: 



TZ~"{q,p) = 2tt(5(p -q >/ + sign(p ) / d 2 x i 



U(x_ 



\sign(p 



(11) 

3 A word of caution is required here. The vacuum-to-vacuum factor is necessary for unitar- 
ity, and this approximation can lead to probabilities that are larger than one. This happens 
because the "smallness" of (Za) 2 can be overcompensated by the large logarithms that ac- 
company it. However, this will not be a problem as long as the multiplicity n remains much 
smaller than unity, as discussed later. 
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with 



U(x±) = Tcxpj-^y dz ^-p a (z ,x ± )t a ^ , 



(12) 

_L 

where the t a, s are the generators of the fundamental representation of SU(N). 
At lowest order in Za, the scattering matrix of the quark in the electromagnetic 
field of the first nucleus simplifies to 

iTW) = 2.%+ - g+ y %x-.x).* > (13) 

KQx - P±) 

where b is the impact parameter in the collision of the two nuclei (we chose the 
origin of coordinates in the transverse plane so that all the b dependence goes 
into the electromagnetic scattering matrix). 

From there, it is easy to work out a compact expression for the amplitude 
to create a qq pair by photoproduction in a peripheral collision. Replacing by 1 
the factor (Oi n |0 ou t)^ at this order in Za, this amplitude is 

M® (p, q) = u(q)T® (q, -p)v(p) . (14) 

For the first two terms in Eq. (^), all the integrals over the ± components of the 
momentum transfer can be done thanks to the delta functions present in the 
individual scattering matrices. For the third term, there are two independent 
transfer momenta. Three of the ± integrals can still be performed with the 
delta functions, and the fourth one can be calculated in the complex plane with 
the theorem of residues. We obtain: 

M? ] (p,q) = 

d 2 kj_ qed ocd u(q)ij + (q l — ft + m)i/-v(p) 



: F (px + qx-k x )F + (fc 



(2tt) 2 ^ ^ ' '"^ + v ~ X) 2p+q~ + (q ± -k ± ) 2 + m 2 

. d 2 k ± qed qcd u(q}J-(ft -tf+m)i/ + v(p) 

+ ( P x + qx-k ± )F_ (kx) 2p _ q+ + {p± _ k±)2+m2 

d2kl ii k ^ F QED ( P x + qx - k 1± - k 2± ) F r(Kx)F« CD (k 2±: 



(15) 



(27t) 2 (2tt) 2 

^(gy+^-ffi + ™y_(fe -rf+m)tJ+v(p) 

X 2q-((p ± - fc 2± ) 2 + m 2 ) + 2p-((qx - K ± ) 2 X 
where we denote 

F2 CD (k ± )= J d 2 x ± e ik ^(u ±1 (x ± )-l} . (16) 

At this stage, it is a simple matter of algebra to combine the three terms into a 
single one: 

M®(p,q) = 
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(PklX d 2 k 2± OED 

I I - ,o - ^ F (p± + qx - «i_ 



(2tt) 2 (2tt) 2 



where we define: 

M(p,q\k 1± ,k 2± ) = 
= u(q) 



U(x ±± )U\x 2± )-l M(p,q\k ±± ,k 



(17) 



2g-((p x -fe 2± )2 



2p~((<Z_L - 



-«(p) (18) 



This is the formula we will use in the following in order to compute the difff active 
cross-section for the photoproduction of a quark-antiquark pair. 



3.2 Equality of Pf cl and n at order (Za) 2 

Before going to the detailed calculation of the diffractive cross-section, it is 
interesting to check a property that holds at order (Za) 2 . At this order, we 
have: 

n = 0{{Zaf) 
pinci = 0((Za) 2 ) 



P*» cl = for n > 2 



(19) 



where P™ cl is the inclusive probability to produce n pairs in a collision. Since by 
definition n = P' ncl + 2P™ cy + ■ ■ •, we immediately conclude that at order (Za) 2 , 
one has n = P{ . Although this property is obvious at such a formal level, it 
is interesting to see how it comes about in the explicit calculation. Indeed, we 
know that n is given by the square of a retarded amplitude, to which only the 
first two diagrams of figure [| contribute |37j] , while all three diagrams contribute 
to P-J ncl . The contribution of the first two diagrams to the retarded amplitude 
was calculated in 13711 , 



Mf (p,q) 



d 2 k± qed 
(2^P 



(PJ 



xu(q){ 



</+(</- 



+ q±- kx) 

m)if-v(p) 



d 2 x<> ik±x± 



i-( 



U(xx) - 1 
tf+m)J + v(p) 



2p+q + (qx — kx) 2 + m 2 



2p-q+ + (px - kx) 2 



r}«(p) 



(20) 



After some algebra, this expression can be rewritten as 



Mf(p,q) = 

. f d 2 k 1± d 2 k 2 x ^Q ED , 



fcu - k 2± ) 
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d x 1 ±d x 2 ±e 



U(x 1± ) - U{x 2± ) M(p,q\k 1± ,k 2± 



(21) 



Therefore, we see that at this order in Za, the time-ordered amplitude ( p"7| ) 
and the retarded amplitude differ only by a factor U(x 2 ±) in the integrand. 
When squaring the amplitudes in order to get either n or Pf cl , this results m 
a factor U(x 2 ±)U^ (x 2 ' ± ). Integrating over the transverse momentum p± of the 
antiquark yields 8{x 2 ^ — x 2± ), and this factor becomes unity (U is a unitary 
matrix) , leading to the expected equality n — p^ ncl . 



3.3 A useful inequality- 
There is another very important property that we can prove without performing 
explicitly all the integrals. Performing only the integrals over p± and q±, one 
can see that the eikonal color matrices present in the square of the time-ordered 
amplitude are of the form 

( [U(x 1± )rf(x a± ) - l] [u(x 2± )U^(x 1± ) - l] )_ 

l-(U(x 1± )tf(x* ± )) z 



(22) 



in the inclusive probability, and 

U(x lX )rf(x aX )-l 



U{x 2± )U\x,_ 



- 1 



l-(U{x xl _)U\x 2i _)) 1 



(23) 



in the diffractive probability. In the above equations, we have used the fact that 
the corr elator (U(x lA _)W(x 2 _ L ))^. is a real number, and symmetric under the 
exchange of the labels 1 and 2. All the other factors in the integrand are the 
same in both probabilities. 

Recalling now the fact that < (U{x- lX )U^{x 2X ))^ < 1 (see the appendix 
of [p7|), we also know that 



l-(U(x l± )rf(x aX ))-. 



< 



i-(u(x 1± )u^x 2± )y p 



From this it is easy to infer the inequality 



1 



pam < _ pi 

i - 2 1 



incl 



(24) 



(25) 



which will provide a useful consistency check at the end of the calculation. Using 
the fact that (Ufaxx)!/' (x 2 x))~ vanishes in the limit Q s — * +oo, we also obtain 
from the above considerations 



pincl 

lim pdiff = £i_ 

Qa->+00 2 



(26) 
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Therefore, the upper bound for Pf lS (the 'black-disc limit') is reached for infinite 
Q s , i.e. in the asymptotic limit where the center of mass energy of the collision 
is infinite. 



4 p± spectrum for diffraction 

We are now in a position to derive an expression for the quantity Pj 3 . For this 
diffractive probability, we need to evaluate the average over the sources p of the 
amplitude M±(p,q). This involves the average (U(x 1± )U^ (x 2± ) — l)~, which 
was calculated in the appendix of p7| , 

(U(x 1± )U\x 2± ) - 1)_ 

= r(x 1± )V(x 3± )(e- B ^ x ^- x ^ - 1 



+(V(x 1± ) + V{x 2l _) - 2P(x 1± )P(x 2± )) (e- Sl - l) . (27) 

Here the function V(x±) describes the transverse profile of the nucleus (see 
[ p7| for more details on how it is introduced in the formalism). We assume 
V(x±) = outside the nucleus, and V(x±) = 1 inside, and that for large nuclei 
the transition at the surface of the nucleus occurs over a distance much smaller 
than the nuclear radius R. In the above equation, we denote: 

B 2 (^±) s Q'i J<fz ± lG„(z ± )-G„(z ± -z: ± )} 2 



47T V^-L^- 



QCD 



Bi = Q' 2 a I d 2 z ± G (z ± ) 2 ~ » 1 (28) 



A 2 

QCD 

with Gq(x±) the 2-dimensional free propagator^ and where Q s is the saturation 
momentum. In the previous equations, one has to introduce by hand the scale 
A Q g^as a regulator for infrared singularities. More details on this can be found 

In Eq. (pTj), the second term depends on x 1 _ j _ and x 2 ^ only via the pro- 
file function V . Performing explicitly the integral over these coordinates, the 
contribution of this term to {^M-p'(p, g)^ is: 

i J ^^^^F QE ° (pj_ + q ± - k 1± - k 2± )M{p 1 q\k 1± ,k 2± ) 
x (v(k 1± )(2n) 2 6(k 2± ) + V(k 2± )(2ir) 2 S(k 1± ) - 2V(k lX )V(k 2 



(29) 



4 This propagator satisfies: d 2 Go{X ± - Z ± )/dZ 2 . = 8(X± - Z±). 
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where the function V(k±) is the Fourier transform of the profile V(x±). This 
function is strongly peaked around k± = 0, with a typical width of order 1/R. 
Therefore, both k 1± and fc 2± are at most of the order of 1/R in this integral. 
Since the factor M(p, q\k 1± , k 2j _) does not contain any scale as small as 1/R 
(this amplitude is controlled by the scale of the quark mass to), it is legitimate 
to replace it by the constant M(p, q\0, 0). By choosing k± = k 1± + fc 2j _ as one 
of the integration variables, this contribution can be rewritten as 

, f CPk QED , 

2iM(p,q\0,0) J j—±F (px + q±-k±) 

(r(k ± )(2n) 2 S(k' ± ) - V(k±)V(k± - k' ± ) 



d 2 k' ± 



(30) 



At this point, it is trivial to see that the function of k± appearing on the second 
line of the previous expression is vanishing if we have P(x±)V(x±) ss V(x±). 
This is indeed the case if the transition of P(x±) from to 1 is sharp, i.e. if 
we can neglect surface effects in the description of the nucleus. To recapitulate, 
one can drop the second term of Eq. (j27j) if wc have to ^ 1/R and R ^ A~* D . 

For the remaining term, it is convenient to resort to an approximation al- 
ready used in Q, valid when Q s 3> A QC£J 3> 1/R: 

T{x 1± )V{x 2± )e- B2{x ^- x ^ ] « V{x 1± )e- B2{x ^- X ^ ) . (31) 

This approximation just indicates that due to the factor exp(— B2) the difference 
x ± j_ —x 2 ^ is much smaller than the nuclear radius. Additionally, using the same 
arguments as above, the term V(x 1 j^)'P(x 2 j^) can be replaced by a single factor 
P(x ± j_). The Fourier transform of the first term of Eq. ( |27| ) is then trivial, as 
we have 



/ 



= V(k 1± + fe 2X )(C(fe 2 x) - (2ir) 2 S(k 2l _)) , (32) 



where, like in j37j, we denote 

C(fcx) = / d 2 x x e ik ^ (U(0)rf(x ± ))-. . (33) 



In the following, we use the notation 

D(k ± ) = C(k±) - (2Tr) 2 5(k ± ) . (34) 

To summarize, we have so far obtained a compact expression for the average 
of M[ p] that reads 

nMt \\ ■ f d 2 k 1± d 2 k 2± qed 

M[\p,q))^ = l J {2it)2 {27t)2 F ( P± + q ± -k 1± -k 2± ) 

xV(k 1± + k 2± )D(k 2± )M(p, q\k 1± ,k a± ) . (35) 
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Note that the factor V(k ± j_ + k 2j _) implies that the total momentum transfer 
k % x + from the nucleus acting via its color field is at most of order 1/R, a 
feature that was to be expected for diffraction. 

From here, we can proceed as in |37|. Since the photoproduction cross- 
section is enhanced by a factor Z 2 for photons that are produced coherently 
by the first nucleus, i.e. for photons with a transverse momentum p± + q± — 
fe 1± — fe 2 ± smaller than 1/R, we can perform a Taylor expansion of the factor 
M(P, q\ki±,k 2± ) around the point where k 1± + k 2± = p± + q±. Up to terms 
of higher order, this gives 

M(p, q\k lX , fe 2X ) « W^M + (p ± + q± _ kix _ feaX ) . g | fe ) ( 36 ) 
/' • '/ 

with 

1 



i(p,g|fc 2 _L) = 7— 



(p~ + g-)(m 2 + (fe a _L - p ± ) 2 ) 

x% _ , _ v+(fc 3 -L -px) + o r« (p) • 

+ q 2 j 

(37) 

Thanks to the sum rule 5 

d 2 k ± 



D(k ± ) = , (39) 



the term of order zero in the expansion ([56|) does not contribute. Integrating 
over the photon momentum m± = p± + q±_ — k 1± — k 2± , while keeping k 2± 
fixed, leads to 

'M®(p, q )\ = 2iZae q V(p x+ q x ) j *^ D(kx) ±I&^ . (40) 
Using then 

]2u b-L(p,q\k ± )b-L*(p,q\k' ± ) _ L{p,q\k ± ) ■ L*(p,q\k' ± ) f d 2 b 

b 2 b 2 2 7 b 2 [ ' 

and 

V 2 {pi_ +q±) ^nR 2 {2ix) 2 8{p 1 _ +q x ) , (42) 

we can write the diffractive cross-section for the photoproduction of a qq pair 
as 

-y/r, 



dm _ 



j d 2 bPf s = nR 2 2N c (Za) 2 e 2 q J ^ 

2R 



5 In order to prove this sum rule, one has to notice that 

d 2 fc, „,, , / t , a r d 2 k 



I W c{k±) = (W®)-, = 1 = / (2 ^ (fe - } ■ (38) 
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d 3 p 



d 3 q 



-(27r) 2 ,5(p ± + <u) 



(2tt) 3 2w p (27r) 3 2w g 
X / |^|^^(^)W)Tr(i(p,9|fe±)-i*(p,g|fcl)).(43) 

In the integration over the impact parameter, the lower bound is determined 
by the geometry of a peripheral collision of two nuclei of radius R, and the 
upper bound is the limit beyond which the electromagnetic field does not have 
a Weizsacker- Williams form. The calculation of the Dirac's trace leads to 



Tr(L(p,q\k ± )-L*(p,q\k' ± )) 



16p q 



(p- + q-) 2 {m 2 + (k± - p ± ) 2 )(m 2 + (k' ± - Pl _) 2 ) 



x |m 2 + (fc ± - Pl _) ■ (k' ± - p±_) 



V 



p~ +q~ 



(p-+q-) 2 



]}■ 

(44) 



At this stage, it is straightforward to perform the integrals over p ± , q and 
q± in order to obtain the cross-section per unit of rapidity (we make use of 
dq+/q+ = 2dy): 



d<7<? 



cliff 



dy 



= ttR 



7/m 

AN c {Za) 2 e 2 r ,n 



r <¥b r d z P 
J ~V J (2^ 

2R 



d 2 P± 



f d 2 k 1 _ d 2 k'^ 
J (2tt) 2 (2ir) 2 



D(k±)D(k' ± ) 



m 



\{k ± -p ± ) ■ (k' ± -p ± ) 



(m 2 + (k ± - p ± ) 2 )(m 2 + (k' ± - Pl _) 2 ) ■ 

(45) 



One can then note that the function D(k±) depends only on the radial variable 
k±, and perform analytically all the angular integrals, which leads to the fol- 
lowing formula for the p± spectrum of the quarks (or, equivalcntly, antiquarks) 
produced diffractively: 



daf s D2 N c (Za) 2 e 2 



■y/r, 



dydpA 



2ir 5 



fb 

'b 2 



2R 



y-pi. 



m 2 ll + - T {h-2p\h 



(46) 
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where we denote^ 



+ OC 



h = I dkj_k ± D(k ± ) 



o 

+00 



dk±k ± D(k ± ) 



\/(m 2 + p\ + k 2 ) 2 - kp\k\ ' 



m + P±_ + fe J_ 

^{m 2 + p'j_ + k 2 ± ) 2 - ip 2 ± k'j_ 



(48) 



Those two integrals are easy to evaluate numerically in order to obtain the p± 

le+03 
le+02 



p 

9 le+01 



1 

-3 
-3 



-3 



le+00 
le-01 
le-02 
le-03 
le-04 
le-05 



0.1 



1 10 

P/ A QCD 



100 



Figure 4: Differential diffractive cross-section as a function of the transverse 
momentum of the quark. The value of Q S /A QCD is set to 10. Open circles: 
strange quarks. Filled squares: charm quark. Open squares: bottom quark. 
The differential cross-section is in units of mbarns/A QCD . 

spectrum of the produced quarks. Results are plotted in figure |]for a value of 
Q S /A QCD — 10 (i.e. Q s w 2GeV, if we assume A QCD sa 200MeV) for strange 
(m/A QCD ps 1), charm (m/A QCD w 8) and bottom {m/A QCD « 23) quarks. 
Corresponding to LHC collision energy, we chose a Lorentz factor of 7 = 3000, 
which is in the regime where the leading logarithm approximation is valid. 



5 Integrated diffractive cross-section 



3 Note that after integrating over the angle, D as a function of the radial variable kj_ reads: 



D(k ± ) = C(k ± ) - ^S(kx) 

K± 



(47) 
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5.1 Cross-section 



The next step is to perform the integration over the transverse momentum 
of the produced quarks, in order to obtain total cross-sections or production 
probabilities. This can be done numerically, and the results are plotted in figure 
Has a function of Q S /A QCD , for strange, charm, and bottom quarks. Again, the 
Lorentz factor is taken to be 7 = 3000. We can see that the cross-section for the 



-3 
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le+02 



le+01 



le+00 



le-01 



le-02 



4> o o 



10 



15 20 

Qs /A QCD 
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Figure 5: Diffractive cross-section in mbarns per unit of rapidity as a function 
of the saturation scale. Open circles: strange quarks. Filled squares: charm 
quarks. Open squares: bottom quarks. 



production of cc and bb pairs is very sensitive to the value of the saturation scale 
Q s , especially in the range Q S /A QCD G [5, 15]. Therefore, this process could be 
used as a way to measure the saturation scale for ultra-relativistic heavy nuclei. 



5.2 Unitarity and black-disc limit 

It is interesting to investigate explicitly the inequality (|2|) between the diffrac- 
tive qq production probability Pf 1 and the inclusive probability P^ ncl (which at 
leading order in Za is equal to n calculated in (37) ) . The probability can easily 
be obtained from the expression for the respective cross-section by integrating 
over the rapidity, but without performing the integration over the impact pa- 
rameter. Note that assuming 7 — ► 00 in our present approach leads to a flat 
rapidity distribution, and a cutoff in rapidity must be introduced by hand. For 
the numerical calculation, we assume^] a rapidity range of 2 111(7), which leads 

7 The precise value of the cutoff does not play any role for the verification of the inequality 
between P^ lfI and Pj ncl /2, it has only an influence on the size of the unitarization corrections. 
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to the relation 



P? iS (2R) = ln 7 



AttR 2 In • 



7 



2mi? 



diff 



(49) 



between the probability at the smallest possible impact parameter and the cross- 
section. In figure ||, we have displayed n, n/2 and Pf lS , at b = 2R and for 
a Lorentz factor 7 = 3000, as functions of Q s / A QCD for strange and charm 
quarks. As expected, we always find P± lS < n/2. Wc also observe that the 



ie+00 



ie-01 



ie-02 



ie+00 



ie-01 



ie-02 




15 20 

Qs /A QCD 



7~1 




le-03 u ' 1 1 ' 1 

5 10 15 20 25 30 

Q S /A QCD 

Figure 6: Verification of the inequality Pi lS < n/2 for the production of qq 
pairs at b = 2R, as a function of Q S /A QCD . Top: strange quarks; bottom: charm 
quarks. Filled squares: average multiplicity for inclusive events (n); solid line: 
n/2. Open squares: non-unitarized probability for the production of a qq pair 
in a diffractive event (P± ); dashed line: naive unitarization of this probability 
(P 1 diff exp(-P 1 diff )). 

two quantities tend to a common value as Q s becomes large. Indeed, in the 
McLerran-Venugopalan model the black-disc limit corresponds to an infinitely 
large value of the saturation scale. 

We have already emphasized the importance of the prefactor (0 ou t|0i n )~, 
to ensure that the obtained probabilities are smaller than 1. However, in the 
present calculation, in which we have kept only the terms of lowest order in Za, 
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this factor is approximated by 1. Therefore, it is expected that this approxima- 
tion violates unitarity, and that unitarity is restored by terms of higher order in 
the electromagnetic coupling. These corrections are rather difficult to calculate 
exactly. For a rough estimate of their size, we suppose that the diffractive mul- 
tiplicities have a Poissonian distribution^, and that the unitarized probability 
to create one pair diffractively is given by P 1 diff exp(-F 1 diff ). This expression is 
also plotted in figure [| We see that the unitarization effects (i.e. the higher 
order corrections in Za) are not important as long as the average multiplicity is 
much smaller than one. In particular, we estimate that the unitarity corrections 
remain small in the range of values of Q s expected to be relevant at RHIC or 
at LHC (between 1 and 3GeV). We also emphasize the fact that this estimate 
has been done for the probabilities at the smallest allowed impact parameter 
b = 2R, where the non-unitarized probabilities are the largest. Integrating the 
unitarized probability over the impact parameter for an estimate of the unita- 
rization corrections in the diffractive cross-section, they are found to be smaller 
by a relative factor of P dlff (2R)/2 \n('j/2mR). Hence, they can be neglected 
here. 

5.3 Summary 

To close this section, let us provide a table of the values of the cross-section 
per unit of rapidity, for charm and bottom quark pairs, and for inclusive^] and 
diffractive production: 





charm 


bottom 


da iuc i/dy (mb) 


355 


11 


da d is/dy (mb) 


60 


0.62 



Those numbers are obtained for a Lorentz factor 7 = 3000, a nuclear radius 
R — 6fm and a saturation scale Q s = 10A QCD . 



6 Conclusions 

In this paper, we have evaluated the cross-section for the diffractive photopro- 
duction of qq pairs in peripheral collisions of heavy nuclei. Diffractive events 
are particularly interesting because they could provide a cleaner final state due 
a rapidity gap between the nuclei and the pair. This process can be used as a 
way to test saturation models, and to measure the value of the saturation scale 
in heavy ion experiments. Indeed, the yield of qq pairs has been found to be 
very sensitive to the value of the saturation scale, especially for c and b quarks 
and for values of the saturation scale between 1 and 3GeV. 



"This supposition do ns n ot hold strictly due to correlations in the amplitude to produce 
more than one pair, see fej] . 

9 The inclusive values come from [p7[ , corrected by the appropriate factors of for the 
electric charges of the quarks, which were not included in [M . 
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We have also checked that our result is consistent with the one obtained 
in |37j for the photoproduction of pairs in inclusive events. In particular, the 
inequality P^ lff < P 1 lncl /2 is satisfied by our calculation, and the black-disc limit 
is recovered when the saturation scale becomes infinite. 
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